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ABSTRACT 


A Boolean function may be expressed, using Exclusive 
OR, AND, and NOT gates, in an unique canonical form, called 
Reed-Muller expansion, which is similar to rainterm expansion 
of a Boolean function. The expansion with Exclusive-OR gates 
possesses inherent advantages over conventional realisation 
of Boolean functions* Minimisation of Exclusive-OR gates, 
to suit optimal design conditions, consists of obtaining what 
is known as the polarity function of the coefficients in Reed- 
Muller expansions; polarity function give^ logical conditions 
of the input variables. Transformations of minterm coeffi- 
cients to the coefficients in Reed-Muller expansions, being 
linear, are represented by matrices, called Reed-Muller 
matrices. Reed-Muller matrices are factorised using direct - 
product Good's techniques, which facilitates handling of 
these matrices in the computer. Algorithmic methods for 
finding the eigen-structure of Reed-Muller matrices, using 
the concepts of the null-space, generalised inverses of 
matrices, have been evolved and : . utilized for the genera- 
tion of coefficients in Reed-Muller expansions from which 
one can obtain polarity functions for the minimum 
expansion. 
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CHAPTER 1 


INTRODUCTION 


1.1 SCOPE OF WORK: 

An attempt has been made in this thesis to study 
the eigenstructure of Reed-Muller matrices and utilize 
their structure for generation of Reed-Muller expansion 
coefficients of Boolean functions. An alternative procedure 
based on factorization of Reed-Muller matrices, using Good’s 
techniques, has also been given for the 
generation of Reed-Muller expansion coefficients. 

Any Boolean function may be expressed in a sum(s)~ 
of-product canonical form. Employing only Exclusive-OR 
(Ex-OR) gate and AND gates and by algebraic substitution and 
Simplification this may be uniquely represented by a similar 
canonical form [l],[2]. As an example, a 3 variable function 
of (x 2 »X 2 »Xj,) is considered. The m interm expansion , 
using AND gates, OR gates and NOT gates of this function is 


f ( X 3 ,X 2 ,Xi)=Cq X 3 X 2 X 1 + ^2 ^3 ^3^2^1 


+ C4 X3X2Xj^+C5 X3X2Xj^+ c^ ^ 3 ^ 2^1 ^7 ^^ 3 ^ 2 ^! 


( 1 . 1 ) 
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where c. fo,l? depending on the function value of the related 
input variables. Every OR gate in the above expression may 
be replaced by an Ex-OR gate as only one minterm becomes 
1 for every combination of values x^^x^jX^ and others 0. The 
equivalent expression, therefore, may be written as 

|,(x3,X2,Xi) = c^ x^x^Xj^ (+) c^ x^x^x^ 0 c^ x^x^x^ (i’ ^3 ^ 3 ^ 2 ^! 
©^4 © *^5 ^ 3 ^ 2^1 ® ^ 3 ^ 2^1 

( 1 . 2 ) 

__ * 

and substituting the identity x^ = 1 0) the expansion 

is obtained in the following form 

|^(x3,X2,x^) = c^(l© X3)(1@X2)(10 x^)@ c^(1@X3) 

(1@X2) Xj^Qc^Cl 0 X3)x2(l©Xj^) 0 03(10x3) 

Vl®‘=4 ^3^^®''2^ ’'l 

(3)c6 ^3^2^^©^!^© ^7 ^3 ^2 ^1 

(1.3) 

which reduces to 

= S@^S®^l^’'l©^S®^2^ '^2 
@^S©''l®''2@‘=3)Vl®^S®''4> ^3 
©(‘^o®‘^l©^4©^5^ ^3^1©) 
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^2 ® ^3 ® ^4 © ^5 ® ^6 ® ^ 7 ^ ^ 3 ^ 2^1 
“ ®o © ®1^1 ® ®2^2 © ®3^2^1 © ®4^3® 

^5H^I @ ^6''2©^7''3^2^1 ^ 


where 

% = S 
= %®^1 
02 - c^ e <=2 

®3 0 ^2 ® ^3 

34 = Cq 0 C4 

®5 “ ^o ® ^1 ^ ^4 ® ^5 
®6 “ <=o©^2©^4G^^6 

^7 = Co®«^l® ''2©''3©''4®''5® ‘'6®‘^7 


(1.5) 

and a. So,lu are constant coefficients that are determined 

^ L J 

by the value of the function for each combination of values 
As there|jio complemented variables in (1.4), 
this is also termed as a polynomial expansion. This expansion 
is unique because the expansion with Ex-OR and minterm is 
unique and the conversion to the expansion with a^s as the 
coefficients is also unique. 
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Thus any m variable Boolean function may be expressed 
as a polynomial expansion. These are also known as positive - 
canonic or complement free Reed-Muller (CFRM) expansions of 
2® terms [3], [4], In general, the expansion may be written as 


f (Xi,X2, © aj_x^ ® a^x^ (+) . . . **’ ’^m 


( 1 . 6 ) 

By complementing any variable a non-positive canonic expansion 
also called as the generalised Reed-Muller (GRM) expansion 
may be obtained, and is written as 


f (x 




” *0® Vl©^2*2© 


® 2 ®- 1 ^ 1 ^ 2 ’^ 3 ' * 


(1.7) 


where each x^ input variable may appear in positive or negative 
form but not in both forms. 

Hence, an m variable Boolean function may possess 2® 
different positive and non-positive canonic RM expansions, 
depending on the polarity of the input variables, which may 
be called as the polarity function of the variables of the 
function. Using the following identities it is' possible to 
obtain any non-positive canonic expansion from either the 
positive canonic or any other non-positive canonic RM 
expansions. 
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= l(+)Xj^ and = l(3^i 

101=0 x ^@ x ^=0 ( 1 . 8 ) 


Suppose the positive canonic expansions in 3 variables is 

; 

given as 


f(x3,X2,x^) = a^©a^x^@a2X2@a3X2Xj^0a^X3 


then the expansion with polarity function 110 may be oDtained 
by the substitution of = 10 x^, and therefore obtain 

f(x3,X2,Xj^) = (aQ@a;L)©a^Xj^@(a2 0a3) X2@a3X2X^ 

®7 ^ 3 ^ 2 ^^! 


( 1 . 10 ) 


Coefficients 

Where a^s are the ^ of the positive canonic RM 

expansions. Allowing the coefficients of this expansion 
to be de-noted as , (the decimal superscript devoting the 
decimal value of the binary, which represents the variables 
being complemented) these coefficients may be expressed in 
terms of the minterm. coefficients as 

= %© 

ai*- = aj^ = c^Q) 
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*2^ “ *2® *3 
= ®3 

* 4 ^ * ® 4®®5 
* 5 ^ = ®5 

~ ®6®®7 

Z3 51^ 


= =1 ® =3 

= =o®''l®'2® “=3 

= ‘=1©'=5 

= %©''1®<=4©''6 
= '= 1 ©= 3 ©‘= 5©'=7 
= '=o®=l©'=2©‘=3®=4 
©= 5 ©'= 5®‘=7 

( 1 * 11 ) 


The constant coefficients a-^s of the positive/non- 
positive canonic RM expansion are respectively called as 
positive/non-positive canonic FiM expansion coefficients or 
RM/GRM expansion coefficients also. The constant coefficients 

c^s, of the minterm expansion is accordingly, called as the 

or 

minterm expansion coefficient [_ merely minterm coefficients. 


The classical problem in switching theory has been 
the realisation of a switching function with minimum number 
of logical circuit elements, Karnaugh map and QuinenMcClusky* s 
algorithms have helped obtain a limited minimum conditions. 
Besides these algorithms, decomposition and other techniques 
have been employed to achieve economical multilevel expressions 
for switching functions [5j, [6]. Geometric formats like 
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the mapping of the positive-canonic RM expansion coefficients 
are used to generate any Ex-OR realisation [7] and thus try 
and obtain a minimum Ex-OR design. Boolean realisations use 
AND/OR/NAND and NOR gates, whereas these same functions can 
also be implemented using Ex-OR/AND/NOJI; gates. 

The advantages of Ex-OR logic design are that more 
economical realisation in terms of the number of gates and 
the number of gate interconnections are possiblej also, designs 
with Ex-OR gates facilitate easier testing in comparison 
with vertex networks since the change of any input to an Ex-OR 
gate will propogate a change to the output, unlike the 
the vertex Boolean realisation which require specific input 
pattern changes to sensitise a path to the output [8]. 

However, the testability advantage may even exceed the 
possible disadvantage of the cost factor in the Ex-OR 
realisation. 

The mapping of the minterm coefficients onto the 
RM expansion coefficients may be considered as the mapping 
of the vector space C to a vector space A over the binary 
field F. As this mapping may be written as a set of n equations 
in n variables viz, the a^s and where n = m being the 
number of variables in the Boolean- function, the mapping may 
be considered as a linear transformation TtC-^A, ^ich 
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satisfies the condition (g c^) T = g( c^T) (J)h( c^T) , 

fox all vectors c^ and c^ in C and all scalars g,h in F, viz. 
0,1. If basis of the vector space C and 

^1*^2* **^n vectors in A, then there is one and only 

linear transformation T:C-:?‘A with Cj^T = a^^, ..^c^T = a^ [ 9 ]. 
Allowing i^ to be the unit vectors, ij^ = (1,0,... O), ... 
in = (0,0, ...1) of the vector space G, then each may be 
written as 

V = = (^11-^12 ^in) 

IjT = 


i T » 
n 


n 


* (i 


nV n2 


» • • 


•^nn^ 


( 1 . 12 ) 


and therefore there is just one linear transformation 
associated with (l.,12) and this transformation is determined 
by the nxn matrix R = (r^^). Therefore there is a one-to-one 
correspondence between the linear transformation T;C-^A and 
the nxn matrix R, with entries from the binary field F. 

Given T, the corresponding matrix R is the matrix with ith 
row, the row of co-ordinates i^Tj given R=(r. .), T is the 
unique linear transformation carrying each unit vector i^ 
of C into the ith row ^in^ 
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Therefore, the processing of the minterm coefficients 
c^s may be expressed in the matrix form; by allowing the 
minterm coefficients and the GRM expansion coefficients to be 
column vectors; this may be written as 

a = Rc (1.13) 

where R is a constant matrix, called as the Reed-Muller 
matrix, with matrix multiplication operation being modulo 2, 
as all operations are over GF(2). The example of the 3 
variable Boolean function may be written as 


‘^0 
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The RM matrix for the positive canonic RM expansion 
may also be iteratively built up from [lO]. 



[ 1 ] 


(1.14) 


where m’ = m~l, for the m variable Boolean function. 

In a similar fashion, it must be feasible to obtain 
the RM matrices for the non-positive canonic RM expansions. 

The scheme developed is not an iterative one as in ( 1.14) 
but the matrix is obtained as a direct-product of two kinds 
of matrices of order 2x2. These matrices are termed as core 
•matrices. The core matrix associated with a positive polarity 
variable is the 2x2 matrix and the one associated with 

the complemented variable is . RM matrices for the 

non— positive canonic RM expansions are thus obtained as a 
direct-product of a combination of these two kinds of core - 
matrices, combined as per the polarity function of the 
variables in the Boolean function. RM matricies may also be 
factorised using Good’s technique. 

Based on the manner in which the RM matrices may be 
factorised two classes of RM matrices may be had viz., the 
RM matrix for the positive canonic RM expansion and for the 
non-positive canonic RM expansions. Salient features of the 
RM matrices observed are: 
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i) the FIM matrix for the positive canonic RM expansion 
is involutary, 

ii) the roots of the characteristic equation of the RM 
matrix, for the positive canonic RM expansion, is 
in the GF(2) 

iii) the roots of the characteristic equation of the PM. 
matrix is repetative 

iv) the characteristic equation of the RM matrix/ for the 
non-positive canonic RM expansions, is a product of 
second degree irreducible polynomial, 

v) the roots of the characteristic equation” of the RM 
matrix, for the non-positive canonic RM expansions, 
are therefore to be found in the extension field 
also, 

vi) the number of kinds of roots being limited to the 
elements of the extension field, the roots of the 
characteristic equation of the RM matrix, for the 
non-positive canonic RM expansion, is also repetitive. 

The characteristic roots are found from the characteri 
Stic equation of the matrices, and therefore, if characteri- 
stic roots are obtained as direct product of eigenvalues of 
the core-matrices, from which the RM matrices are constructed 



12 


then it is possible to construct the characteristic equation 
from these characteristic roots. The characteristic roots are 
known as eigen-values also. 

Every eigen-value has an eigen-vector corresponding 
to it. When the roots of the characteristic equations are 
repetitive, unique eigen-vectors for each repetitive root 
can not be found by the normal procedures. Two methods of 
obtaining the modal matrix, the matrix whose columns are the 
set of linearly independent (LIN) eigen-vectors of a matrix, 
have been outlined. In one, the concept of generalised 
inverses of matrices is used, and in the other the concept 
of null-space of a matrix is employed. The second method has 
been suitably adopted, for the RM matrix for the positive canonic 
RM expansion, and a simplified procedure has been evolved. 

Having obtained the eigen structure of the RM matrices, 
algorithms to generate the expansion coefficients from the 
given truth-table of a Boolean function, have been evolved. 

The algorithm for the positive-canonic RM expansion makes use 
of the eigen-vectors of the RM matrix and relates all 
possible minterm coefficients to its corresponding expansion 
coefficients. The set of eigen-vectors forms a subgroup of 
the group of all possible minterm coefficients and therefore, 
cosets may be formed and these minterm coefficients may 
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now be associated to the coriesponding expansion coefficients 
through the eigen-vectors. The other algorithm makes use of 
the eigen-values for generation of the non-positive canonic 
RM expansion coefficients, Good’s technique for factori- 
sation of RM matrices may also be used for the generation of 
the RM expansion coefficients, 

A geometric mapping of the positive canonic RM 
expansion coefficients, called as ’b. maps’ , has been 

vl 

developed to generate all non-positive canonic expansion 
coefficients [?]. However, the map operations become 
difficult when the number of variables in the Boolean function 
increases beyond 6. A computer method has been developed, 
to generate all non-positive canonic RM expansion coefficients 
from the positive canonic RM expansion coefficients [lO]. 

1.2 LINE OF APPROACH: 

The RM matrices wore defined in the previous section ^he 
RM matrix are constructed as direct product of core-matrices. 
Another method of factorising the RM matrix based on 
Good's technique is also considered. They are classified 
as RM matrix, for the positive canonic RM expansion and the 
RM matrices for the non-positive canonic RM expansions. The 
generalisation is based on polarity function of the variables 
in the Boolean function. 
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The eigen-structure of the RM matrices would be 
different for different structure of the RM matrices. It 
has been possible to define the number of eigen-values of 
the RM matrices as a function of the total number of variables 
of the Boolean function and the number of variables of these 
that appear complemented. Hence, eventhough the eigen-values 
of two different polarity functions may be same, as the 
number of variables and the number of variables being comple- 
mented being the same for two different polarity functions, 
the eigen— vectors of the RM matrices with these two polarity 
function would not be the same. The RM matrices are 
obtained, as direct product of the core matrices, in sequence 
as per the pattern of the polarity function and hence the 
RM matrices for different polarity functions would be 
different and so would the eigen-vectors be, 

A Boolean function is expressed in binary and if it is 
not possible to obtain the eigen-values in the binary field 
GF(2), then the roots are found in its extension field. The 
characteristic equations of the RM matrices, for the non- 
positive canonic RM expansion, being a product of an irre- 
ducible polynomial of degree 2, it is only natural to expect 
the roots of the characteristic equations to be repetitive, 
being confined to the eluents of the field. The method of 
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obtaining the eigen-vectors when repetitive roots occur is 
based on the generalised inverse of the characteristic matrices, 
A specific technique has oeen developed in this thesis to find 
the modal matrix of the RM matrix for the positxve-canonic 
RM expansion. 

Eigen-structure of any matrix characterises the 
transform. The eigen-vectors of the RM matrix for the positive 
canonic expansion, partitions the vector space C, the 
vector space of all possible sets of the minterm coefficients 
c^s. Vifith the set of eigen-vector as a subgroup, cosets 
have been formed of this vector space C. By identifying the 
coset-leader and associating the coset-leader with an 
eigen-vector, through the characteristic matrix, thQ minterm 
coefficient is related to its corresponding positive canonic 
RM expansion coefficients. Similarly, the eigen-values of 
the RM matrix, for the non-positive canonic expansion, 
have been associated with the characteristic matrix and the 
minterm coefficients, to generate the non-positive canonic RM 
expansion coefficients. 

1.3 OUTLINE OF Tlffi..THESlS: 

In this chapter the aim of this thesis and the common 
terms used in this thesis were explained. The line of 
approach towards the generation of the RM expansion coefficients 
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was also given* The study of the eigen-structure itself was 
considered important because some input vectors remain 
unaltered by the transformation, except by a scaling factor, 
called the eigen-values. These input vectors are called as 
the eigen-vectors and they characterise the transformation 
and these could be put to appropriate use. 

Chapter 2 is a mathematical recaptulation of the pre- 
liminaries necessary in the development of this thesis and 
highlights some techniques for obtaining the eigen-vectors. 

In the next chapter, Chapter 3, the structure of the RM 
matrices have been studied and has been represented as direct 
product of core matrices. Properties of direct product of 
matrices have been applied to find the eigen-values of the 
RM matrices. The direct-product method of finding the eigen- 
values of the RM matrices, for tho non-positive canonic RM 
expansion, from the eigen-values of the core matrices have 
been simplified into formulae. Characteristic equations of 
the RM matrices, from the eigen-values so obtained, has also 
been outlined. A simple method of finding the modal matrix 
of the RM matrix, for the positive-canonic RM expansion, has 
also been developed in this chapter. To obtain the modal 
matrix of the RM matrices, for the non-positive canonic RM 
expansion, the generalised inverse of the characteristic. 
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matrix, especially when the roots of the characteristic 
equation are repetitive, has been suitably adopted^ The 
FIM matrix has also been factorised using Good’s technique. 

In Chapter 4, the minimisation of the Ex-OR design 
has been discussed. Algorithm for generation of the 
positive canonic RM expansion coefficients and non-positive 
canonic RM expansion coefficients, directly from the minterm 
coefficients have been given. These algorithm are based on 
the eigen-vectors of the BM. matrices and the eigen-valuos 
of the RM matrices respectively, Basea on Good’s technique 
of matrix factorisation sequential (binary/decimal) genera- 
tion of all RM expansion coefficients has also been done. 

The last chapter summarises the work in this thesis. 
Suggestions for further work in this field has also been 
given in this chapter* 



CHAPTER 2 


MATHEMATICAL PRELIMINARIES 

This chapter briefly outlines the mathematical pre- 
requisites, These include relevant results from the theory 
of groups, rings, fields and matrices. Matrices not only 
simplify the analysis but help in organising computer methods 
for the techniques developed. Here the algebra of matrices 
essential for the thesis has also been included. Standard 
results have been adopted from a number of text books, a 
few of these are [9], [12], [l3] , [l4] , [l5j , [l6], [l7] , [l8], 

[19], [20]. 

2.1 GROUPS, RINGS AND FIELDS: 

2.1.1 Groups: 

A group is a set of elements with one operation and 
its inverse, such as addition and its inverse viz. subtraction | 

and multiplication and its inverse viz, division. A group G | 

thus is a set of objects for which an operation is defined j 

and for which the following rules hold. Let a,b, c be | 

elements of the group. The operation is a single valued 

t 

function of two variables, sayf(a,b) = c, which is [ 

denoted as a + b = c or ab = c and called as addition or 
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"multiplication even though it may not oe addition ox 
multiplication of the arithmetic of ordinary numbers. 

Rule l (Closure); The operation can be applied to any two 
group elements to give a third group element. 

Rule 2 (Associative Law): For any three elements a,b and c 
of the group, (a+b)+c = a4'(b+c) if the operation is written 
as addition , or a (bc)=(ab)c if the operation is written 
as multiplication* 

Rule 3 : There is an identity element. If the operation 
is addition, the identity elemat is called zero, 'O’ and 
if the operation is multiplication, the identity element 
is one, ' 1' , 

Rule 4 t Every element of the group has an inverse element. 

If the operation is addition, the inverse element corres- 
ponding to 'a‘ IS '-a' and is defined by the equation as 
a+(-a) = (-a)+a=0. If the operation is multiplication, the 
inverse of ’a' is ' a*”^' and is defined as aa”^ = a~^a = 1. 

In addition to the above laws, a group may satisfy the 
commutative law; i.e. a+b = b+a, or if the operation is 
multiplication, ab = ba. Such a group is called an Abelian 
or a Commutative group. 

2.1,2 Rings; 

A ring R is a set of elements for which tv/o operations 
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are defined. One is addition and the other multiplication, 

even though these may not be ordinary addition or multipli- 

for 

cation of numbers. In order l_ R to b e a ring, the following 
rules must apply: 

Rule 1 ; The set R is an abelian group under addition. 

Rule 2 ; (closure). For any two elements a and b of R, the 
product ab is defined and is an element of R. 

Rule 3 (Associative Law): For any three elements a,b, and c 
of R, a(bc)=(ab)c. 

Rule 4 (Distributive Law); For any three elements a,b and c 
of R, a(b+c)=ab+ac and ( b-tc ) a=ba+ca . 

A ring is called commutative if its multiplication operation 
is commutative, i.e. for any two elements a and b ab=ba* 

2.1.3 Fields: A field is a commutative ring with an unit 
element (multiplicative identity) in which every non-zero 
element has a multiplicative inverse. 

2.3 SUBGROUPS AND GOSETS: 

A subset of elements of a group G is called a subgroup 
H, if it satisfies all the rules for a group itself. H is a 
subgroup if it satisfies the rules for closure and the rule 
of inverse also applies*. For a set which fulfils the two 
rules, the identity must also be present. 
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Suppose the elements of a group G are g2fSf2»^3* * • » 
and the elements of a subgroup H are hj^, h 2 ,h 2 . . and consider 
the array formed as follows: 

The first row is the subgroup, with the identity at 
the left and each other element appearing once and only once. 
The first element in the second row is any element not 
appearing in the first row, and the rest of the elements are 
obtained by operating on each subgroup element by this first 
element on the left. Similarly a third, fourth and fifth 


rows are form 

/ 

,ad each 

with a 

previously unused group element 

in the first 

column, 

until 

all the group 

elements appear 

Somewhere in 

•’i 

the array. 

= Ir ^2 

» ^ 3 * * * * 

- % 

9l \ 


^1^2 

> • • * 

' Sir'll 

92*’i 

• 

• 

= 02* 

^ 2^2 

t ^ 2 ^ 3 * * * * 

’ 32''n 

• 

^m 1 

= 


, g^h3, ... 

' 9m*’n 


The set of elements in a row of this array is called 
a left coset and the element appearing in the first column 
i*'S called the coset leader. Right cosets could be similarly 
formed. Also every element of the group G is in one and 
only one coset of a subgroup H, 
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2.3 GALOIS FIELDS: 

Before defining the Galois fields, associated termi- 
nologies are explained, 

2.3.1 Ideals, Residue Class and Residue Class Rings: 

An ideal I is a subset of elements of a ring R with 
the properties that I is a subgroup of the additive group of 
R and for any element ‘a’ of I and any element 'r* of R> ar 
and ra are in I. Since an ideal is a subgroup, cosets can be 
formed. In this case the cosets are called residue classes 
and all properties of cosets apply to residue classes also. 
The residue classes of a ring with respect to an ideal forms 
a ring and this ring is called the residue class ring. 

2.3.2 Ground Fields, Extension Fields: 

The residue class ring modulo m is a field if and only 
if m is a prime number. These fields are called prime fields, 
or Galois fields of p elements, GF(p). 

If r(x), s(x) and t(x) are polynomials and r(x) s(x) = 
t(x), then it is said that t(x) is divisible by r(x) or that 
r(x) divides t(x), and that r(x) is a factor of t(x). A 
polynomial p(x) of degree n which is not divisible by any 
polynomial of degrees less than n but greater than 0 is 
called irreducible. 
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Let p(x) be the polynomial with coefficients in a 
field F. If p(x) is irreducible in F, that is, if p(x) has 
no factors with coefficients in F, then the algebra of 
polynomials over F modulo p(x) is a field. The field formed 
by taking polynomials over a field F modulo an irreducible 
polynomial p(x) of degree k is called an extension field of 
degree k over F. The original field F is called the ground 
field, 

2.4 GENERALISED INVERSES OF MATRICES: 

The concept of division does not exist in matrix 
algebra, however, an equivalent operation results when the 
matrix is multiplied by an unique matrix, called an 'inverse' 
to obtain the identity matrix, the 'one' of matrix algebra. 
The idea of inverse is that if Ax = b then x = A^^b and this 
A is such that A A = I, This A exists only if matrix 
A is non-singular and is a square matrix. 

This unique inverse is at times denoted as M also, 
as this matrix M satisfies the four conditions given by 
Moore and Penrose and hence at times it is also known as the 
Moore-Penrose inverse. The four conditions are, 

i) ML\ = A 

ii) MAM = M 

iii) AM is symmetric and 

iv) MA is symmetric. 
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The matrix M defined by the four Moore-Penrose 
conditions is unique for a given A. But there are many 
matrices G which satisfy just the first Moore— Penrose 
condition* 

AGA ^ A 

Any matrix G that satisfies this condition is called a 
generalised inverse of A, fox A of dimension pxq and G of 
dimension qxp. Thus many matrices G may be found that 
satisfy this condition, the one exception is when A is non- 
singular in which case there is only one G and it is the 
regular inverse G = A“^=M. The matrix G has numerous other 
names, depending on the various number of conditions it 
satisfies. The M inverse of A is denoted by A*~'^ and the 
generalised inverse by A*". 

2.4.1 An Algorithm to Find the Generalised Inverse G; 

Matrix A has dimension pxq and A^^j^ is the leading 
sub-matrix and is non-singular of rank r^. Then a generalised 
inverse may be obtained as 




"^12 


f-i 

Vi 

0 1 

X q “ 

i 

1 

1 

^21 


is G^ = 

qxp 


£j 


where the null matrices in G have appropriate order to make 
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G q X p, the verification of AGA = a involves using 
^22 ^21 "^11 ^ 12 * 

However if is not a non-singular sub-matrix of 

order r, then by elementary operations on the matrix a 
non-singular submatrix of order r may be brought to the 
leading position and then by partitioning the matrix B, 
and applying the same elementary operation on the generalised 
inverse so obtained, the generalised inverse G of matrix A 
may be obtained. Let R and S be the permutation matrices 
such that RAS brings a non-singular matrix of rank r^ to 
the leading position. Thus 



then F, the generalised inverse of B is F = 

j 

and therefore G = SFR, which is a generalised inverse of A. 

2.5 DIRECT PRODUCT OF MATRICES: 

The direct product of two matrices Ap q and ^ n 
may be defined as 

^pxq® ®raxn 
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and is called as Kronecker product or Zehfuss product. 

The product matrix is partitioned into as many submatrices 
as there are elements of A, each subraatrix being B multi- 
plied by an element of A. Therefore the elements of the 
direct product consist of all possible products of •' elements 
of A multiplied by an element of B. It has, in general, 
order pm x qn. 

Direct product, therefore, is an unique manner of 
factorising a matrix defining a transformation and hence 
implementation of the transformation can be accomplished 
with reduced number of arithmetic operations. 

Some useful results of direct product of matrices are 


a) "hh 

where R - A(0B(^C ... and mixed radices. 

b) If a^^ and are the eigen-values of A a matrix and 

bii and and so on of the others then the eigen-values 
of the resultant matrix a(^b(^.,. is the product of their 
individual eigen-values. 

i.e. ^21^11' ^lAi' ^21*^11 ^21^21 


c) If 

®21 *’21 


are simple roots of the characteristic 


equation of the matrices A and B, then the eigen-vectors of 
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of C = A B are 

^ 11^21 
®2 1*^11 
® 2 1^21 

2.6 EIGENVALUES AND EIGENVECTORS OF MATRICES; 

For u a vector and y a scalar, the equation 
Au = yu 

may be valid for a given matrix A. The conditions for which 
this equation is valid, other than fox u = 0, is that A-yl 
IS singular i.e. jA-yll = 0. The above condition may be 
rewritten as 

(A - yl) u = 0 

A non-null solution for u may be obtained as 
u = [(A - yl)" (A-yl)-l]z 

using a generalised inverse of (A-yl) and where z is arbitrary. 

The condition [A-yll = 0 is called the characteristic 
equation of A. For A of order n, the characteristic equation 
is a polynomial equation in y of order n, with n roots 

i 

denoted by Yj[i Y2» • • some of which may be the same and i 

some of which may be zero. These roots are called latent-roots, | 
characteristic-roots, proper roots, eigenvalues or y-roots. ; 

Corresponding to each y. is a vector u. which satisfies 

JL J* 
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A Ui = y±^± and these vectors are 

correspondingly called latent-vectors, characteristic-vector 
proper-vector or eigenvector. And the matrix [A-ylj is called 
the characteristic matrix. 

The set of all vectors u, which form the columns of a 
matrix is called as a model matrix. Model matrix may also 
be found using Jord^x-A canonical matrix [l2], 

2,6,1 Calculating Eigenvectors: 

Supposing Yj^ is an eigenvalue of A. It is a solution 
of the characteristic equation [a-yII = 0, Determining an 
eigenvector corresponding to Yj^ implies finding a non-null u 
to satisfy Au = which is equivalent to solving 

(A - y^l)u = 0 

The number of linearly independent solution to the above 
equation are n-r(A-Yj^I) i.e, the order of the matrix minus 
the rank of the characteristic matrix (corresponding to the 
characteristic value) and these solutions are 

= [(A-Yj^I)*" (A-Y^D-IJz 

since Yj^ is such that |a-Yj^I| = Oi is singular and so 

at least one non-null solution always exists. For each Y{j» 
the characteristic matrix [a-Yj^I] is processed by elementary 
operations and a matrix is obtained to enable partition 
the matrix as, 



^ ;7 



such that the submatrix Rj^ is non-singular with the same 
rank as the characteristic matrix [A-Yj^I] and then obtain 
the generalised inverse of the characteristic matrix as 



and so obtain the solution 



A number of eigen-vectors Uj^ are found from the 
general solution for each eigenvalue Yj^ of the matrix 
including if if occurs. Since Yj^ o solution to 

the characteristic equation it may be a solution more than 
once, in which case it is called a multiple root. If is 
a solution only once it is called as a simple root. When 
Yj^ is a simple root the rank of the characteristic matrix 
r(A-Yki) n-1. Therefore there is just only one LIN eigen- 
vector corresponding to Yj^* This is reflected in w as w 
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becomes a scalar. Vi/hen Yj^ is a solution more than once 
it is called as multiple roots ond the number of times that 
it is a solution is called its multiplicity. For each 
multiple eigenvalue Yjj» with multiplicity m^^, the rank of 
A-Yj^I must be ascertained as this is the value that plays 
an important role in the determination of oigen~vector s . 


2.7 MJLL-SPACE AND DIMENSION OF NULL-SPaCE OF MATRICES: 

In Sec, 1.1, the linear transformation was identified 
with an n X n matrix R. T was the linear transformation of 
a vector space C into a vector space a. The null-space of a 
linear transformation T is the set of all vectors such that 
= 0, and the null-space of a matrix R is the set of all 
row matrices X which satisfy the homogenous linear equation 
XR = 0. Alternatively, the null-space of the row space of a 
matrix is called the null-space of the matrix. A vector is 
in the null-space of a matrix if it is orthogonal to each row 
of the matrix. The dimension of the null-space of a given 
matrix is called as nullity of the matrix. The rank and 
nullity of an n x n matrix is related by th^ following 
equation. 

Order of the matrix - rank of thematrix = nullity. 



CHAPTER 3 


EIGEN-STRUCTURE OF REED-MU LLER i\/IATRICES 

Methods to obtain the positive canonic RM expansion 
coefficients from the minterm coefficients (obtained from the 
truth-table entries of a Boolean function) have been 
developed [?], [lO], [2l], [22], [23j. In the methods developed 
in [7],[l0] other non-positive canonic RM expansion coefficients 
have been obtained from the positive canonic RivS expansion 
coefficients. In the geometric mapping method [7], the map 
folding operations become difficult when the number of 
variables in the Boolean function increases. Moreover, if 
this difficulty is surmounted then in any map method there 
is no freedom to represent the variables than the available 
3 dimensions, with utmost four variables represented in each 
dimension. The limitations of the b^. maps have been overcome 
in the computer method [lO]. 

In this chapter the RM matrix structures, for different 
polarity functions have been obtained as a direct-product of 
two kinds of core matrices. Techniques have also been 
developed, or suitably adopted, to determine the eigen-structure,. 

3.1 RM MATRICES: 

, c 1 

In Chapter 1 the method of obtaining any non-positive • 
canonic RM expansion from any other expansions by algebraic 
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substitutions and simplification [24] v;as illustratea by 
(1.8), (1.9), (1.10) and (l.ll). For this example the 
transformation, in the matrix form v/ould be 

01000000 
11000000 
01010000 
11110000 
^ ~ 01000100 
11001100 
01010101 
11111111 


This is the transformation from the minterm coefficients, 
obtained from the truth-table entries of the Boolean function, 
to the non-positive canonic RM expansion coefficients. This 
particular transformation pertains to the polarity function ' 

110 of the 3 variable Boolean function. The RM matrix, for 

I 

the positive canonic RM expansion, was also obtained iterativelyi 

I 

(1,14) which may also be obtained as direct product of, what | 

! 

has been termed as the core matrix, for m’ = 0 in (1.14), i.e. i 
the 2x2 matrix [| 9] . To obtain the Rjvl matrix, for the | 

^ ^ f 

non-positive canonic RM expansions, suitable core matrix, [ 

( 

to be called as core matrix for a complemented variable, | 

would have to be appropriately multiplied with the other core j 


I 
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matrices; by multiplication is meant the direct product 
as per some rules. 

The eigen-values and eigen-vectors of the Reed-iVuller 
matrices help in obtaining the positive canonic and non- 
positive canonic expansion coefficients and therefore, the 
eigen-structure of the RM matrices have been developed in 
this chapter subsequently. The results obtained in this 
chapter have been used in the next chapter, in which tne 
algorithms to obtain the expansion coefficients have been 
developed. As the Reed-Muller matrices have been developed 
as direct product of the core matrices, the oigen-values 
too have been obtained as direct product of the eigen-values 
of the core matrices. As brought out earlier the eigen-vectors 
for repetitive roots (eigen-values) of the characteristic 
matrix of RM matrices, have to be uniquely obtained by 

I 

employing the generalised inverses of the characteristic i 

matrix or using the concept of null-space of the- characteristic i 
matrix, | 

3.2 THE CORE MATRICES: ' 

f 

In the previous section the core matrix for the | 

’positive* variable was identified. Like the RM matrix, 
for the positive canonic RM expansion, this core matrix is also i 
involutary, i.e. the matrix is its own inverse. Similarly f 
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to obtain the RM matrices, for non-positive canonic RM 
expansions, a core matrix to represent the complemented 
variable is identified as that matrix which results when 
the core matrix for a positive variable is multipled by the 
counter identity matrix J, This is represented as 


= 


r 

‘•1 1 -* 


and may be called as the core matrix for a complemented 
variable* 

The characteristic equations of the two core matrices 
are respectively. 


(Y + l)^ = 0 (3,2*2) 

(Y^+Y'^'I) “ 0 (3*2,3) 


The roots of the characteristic equation of the core matrix 
is y-1 and y=l. The repetition of a root is also termed as 
its multiplicity and is denoted by mj^. Therefore the 
characteristic roots of the core matrix may be written as 
Y=l, The characteristic equation of the core matrix 

for a Complemented variable is an irreducible polynomial of 
degree 2, Therefore, the roots are not in the ground field, 
GF(2) and are to be found in the first extension field viz. 
(GF2^) and the roots would be from the elements of the 
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extension field, viz. 0,1, a and l+a, with a^-fa4*l = 0 being 
the irreducible polynomial. Hence, the ch: racteristic roots 
of the core matrix of a complemented variable are Y=a and 
Y“l+ct, with multiplicity one each. It is also observed that 
the irreducible polynomial is of degree 2 and therefore the 
characteristic roots always occur as a set. 

3.3 OBTAININ3 RM MATRICES AS DIRECT-PRODUCT OF CORE MATRICES: 


A given m variable Boolean function may be associated 
to the corresponding RM expansion coefficients through trans- 
forms, called as the RM matrices. It was also described in 
Chapter 1, that the RM matrices may be classified as the RM 
matrix for the positive canonic RM expansion or the CFFIM 
expansion and as the RM matrices that represent the other 
non— positive canonic RM expansions, called as the generalised 
Reed-Muller expansions. In this section, the method for 
obtaining the two kinds of Reed-Muller matrices, as direct- 
product of the two kinds of core matrices, is developed. 

3.3.1 RM Matrix for Positive Canonic RM Expansion: 

The RM matrix for the positive-canonic RM expansion, 
in which all the m variables of the Boolean function are 
in the positive form, is obtained as direct product of the 
. The m variable Boolean function is said to 


core matrices 
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have ail positive polarity function and hence each of the 
m variable of the Boolean function is represented by a core 
matrix. Therefore, the Reed-Muller matrix, so obtained 

is of order n, n=2 , For the Boolean function of m variables, 
this may be written as 

^"^^11^ (?) 1 ^ ^ •••• m times (3.3.1) 

The RM matrix so obtained is also involutary. 

3,3.2 RM Matrix for Non-positive Canonic RM Expansions; 

In the non-positive canonic expansion, one or more 
variables of the Boolean function apoear in the complemented 
form, but the same variable does not appear in both the forms 
in the expansion. The RM matrix corresponding to these 
expansions would, therefore, be dependent on the polarity 
function of the Boolean function and is obtained as direct- 
product of both kinds of core matrices. Each variable in the 
function which has positive polarity is represented by the 
core matrix and the variable that is complemented by the core 
matrix of a complemented variable. For the Boolean function 
of m variables, which has the polarity function 1101001... 1, 
the RM matrix is obtained as direct product of the appropriate 
core matrices in that order. 
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R = ® C 


1 0 
1 1 


]fiC° 1 


] ® [ 


1 0 
1 1 


] ® [J 


1 

1 


] ®) [J j;] ® [J: J 


] 


, 0 [ m times (3,3,2) 


Therefore, if there are t variables complemented in the m 
variable functions, there are (m-t) core matrices and t 
core matrices, for a complemented variable. The RM matrices 
so obtained is unlike the RNi matrix for the positive canonic 
RM expansions, i,e, is not involutary. Here too the order of 
the RM matrix is n, 0=2”^. 

3.4 GOOD'S FACTORISATION OF RM MATRICES-. 

The RM matrices have a high degree of redundency and 
if the redundency in the definition of the matrix trans- 
formation is eliminated by matrix factorisation then a more 
efficient means of implementation is possible. Such a 
technique was employed in the development of the fast Fourier 
transform (FFT) and fast Kadamard transform (FHT) and to a 
large class of fast transformations [l4], 

A class of matrices formed by Kronecker product 
operation may be considered as the square submatrices of 
order p, with entries m^ . where i and j range from 0 to 

A » if j 

p-1, and the first index r representing the class of entries 
corresponding to a particular dimension in the Kronecker 
product operation. In general 
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m 


r,o,o r,o,l 
r,l,o 


= 


and 


m 


m_ ^ in_ I 1 . . . ra, 


x,o,p-l 


'^r p-i o ^r, p-1,1* * P~lf P”1 


»1 =“q 


H = 


(3.4.1) 


(3.4.2) 


Hq = ”q-l 


Thus 


H, 


m. 


H, 


loo ^q-1 **’ '^q-l,o,p— 1 M“1 

j"“X j o t q 


H, 


^q-1 “* '^q-l,l»P-l 


m 


q— 1, p**l> o^q“l 


ni, 


H, 


q-l f p— 1 f P“1 


(3.4.3) 


where 


H is a p^ 0 P*^ matrix 
Q 


When operating with Kronecker matrices within a 
computer, it becomes 'desirable to store a representation 
(algorithm) of the entries of the product matrix rather than 
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the matrix itself. The location in the matrix may be described 
by their lexicographic sequence representation. A given 
index of matrix may thus be represented by q digits each 
of which may take the value 0 to p— 1. Representing the 
horizontal index by v and the vertical index by x, the names 
of the rows and columns in lexicographic sequence for H 2 
matrix for p=!3 is 


V 


00 01 02 10 11 12 20 21 22 


00 

X 

01 
02 
10 


H2= 

12 

20 

21 

22 


1 


H2(x,v) 


(3.4.4) 


Representing the v and x variables in the lexicographic number 
system, modulo p^requires q digits to allow v and x to range 


over 0 to 

Therefore v and x 

may be described 

by 

^ = ^n-1 

''n-2 • • • Vo ''1 

p — 1 

(3 

^ ^ ^n-l 

\-2 ••• ’'i 

0 y 1 5 • ♦ • , P — 1 



( 
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Using such notation allows the entries of the p x p matrix 
(3.4,2) to be described by the equation 


H^(x,v) 


p-1 

■n; 

i=o 


p-1 

■It 

j=o 


(x^-i)d(v^-j) 


(3.4.6) 


where d(a-b) is the delta function which takes on th^ value 
1 whenever a=b and 0 otherwise. The representation of 
(3,4,6) may be interpreted as multiplying all entries of the 
core matrix, (3,4.2), together and noting that all but 
one entry would be raised to the 0th power. The- entries of 
the p by p matrix, H 2 » (3,4.4) may now be represented as 


^ 2 ( X , v) 5= 


p-1 p-1 
Tt H 

i=o j=o 


m 


d(xj^-i)(v 

l,i» j 


1 


j) 


P-1 

7C 

i=o 


p-1 

71 

J=0 


d(x -i)a(v -j) 

“o,i,j 

(3.4.7) 


where the exponents determine the correct product of entries 
for a given v and x. In general, the entries for may be 
represented as. 


Hq(x,v) 


q-1 p-1 p-1 d(x^-i)d(v^-j) 

71 11 11 m . 

r=o i=o j=o 


(3.4,8) 


following the recursive notation of (3.4,6) and (3,4,7). 
Representation of the rows or columns of a Kronecker matrix 
in the form of (3.4.8) allows the generation of any single 
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element, column or row of the matrix without storage of the 
entire matrix array. The operation indicated in (3.4.8) 
may be visualised as multiplying all entries of the sub- 
matrices, M^, forming the Kronecker Product, together and 
letting the exponent operation allow only the proper entries to 
be raised to power 1 while all others are raised to the power 

0* For the Kronecker product of identical matrices 

n-1 

p-l p-i I 6(Xj.-l)d(v^-j) 

Hq(x,v) ~ % Ti ? (3.4.9) 

A major observation was made about the matrix facto- 
risation i.e. if highly redundant matrices could be factored 
into a product of matrices with few non-zero entries, then a 
fewer number of operations would be necessary for transfor- 
mation implementation. A technique, attributable to Good 
[26], of matrix factorisation^may be used to decompose the 

class of Kronecker matrices described by (3.4,8) or (3,4.9). 
the 

Thus for Hq Kronecker matrix of (3.4.8), there exists q 

matrices, each of dimension p*^, such that when multiplied 

together will equal H • These matrices may be described as 

9 
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In this matrix there are p*^"**^ non-zero entries and only 
non-redundant elements. Then 




(3.4.11) 


and for the Krocker 


Hn = 


product of identical matrices 
[G]^ (3.4.12) 


The RM matrices obtained earlier (3.3.1) and (3.3.2) 
may be factorised respectively as, 
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where. 


R 


10 


10 


10 


10 


11 


11 


11 


11 


•11 


n columns 


n/2 rows 


n/2 rows 


and, 


R = 


( 3 . 4 . 14 *; 


— _ 


— 

^ 1 

10 


01 


10 

11 


11 


11 

10 


01 


10 

11 


11 


11 

10 


01 


10 

11 


11 

« 


11 

A 

« 

• 


I • 

• 


• 

• 

10 

1 

i 

01 


10 

11 


11 


11 


L 




, as per 
the 

polarity 

function 


( 3 . 4 . 15 ) 


“ ^ times 



45 


where Ra == R and R. = R- as per the polarity function, and 

X. XX 

01 

TJl 

01 

01 

R “ — 

^ 11 

11 

11 

♦ 

11 

both R and R- obtained by rearranging the rows - e€ (3.4.10). 

ilC X 

3.5 EIGENVALUES OF THE REED-MULLER MATRICES*. 

The eigenvalues or the characteristic roots of a 
matrix are the solutions to the equation, 

Ru = (3.5.1) 

for R a matrix and u, a vector and y? ^ scalar. There are 
a number ^of scalars and vectors that satisfies this equation; 
the scalars are called as eigenvalues and the set of eigen- 
values is are termed as the spectrum of R. The eigenvalues 
are also obtained by solving thecharacteristic equation 

( R - Yl) = 0 (3.5.2) 

In the previous section the RM matrices were obtained as direct 
product of the core matrices and were classified as RM matrix 


n/2 rows 

(3.4.16) 

n/2 rows 
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for the positive canonic RM expansion and RM matrices for 
the non— positive canonic RM expansion. The polarity 
function of the variables of the Boolean function determines 
the kinds of roots and their multiplicity. In this section 
the eigenvalues are not found as solutions to the chara- 
cteristic matrices of the RM matrices (3,5.2) but the property 
of direct-product is exploited, 

3.5.1 Eigen-value of the RM Matrix, for Positive Canonic 
RM expansion: 

The RM matrix for the CFRM expansion was obtained as 
direct product of the core matrix (3.3.1); the order of 

the matrix being n=2™. Therefore there are utmost 2® 
eigenvalues. The eigenvalues of the core matrix are 
with multiplicity mj^=:2 and therefore, the eigenvalues of the 
RM matrix for the CFRM expansion are y= 1» mjj^=n. And, 
therefore may be written as 

Spectrum of R = (3,5.3) 

3.5.2 Eigen-values of RM trices, for Non-positive 
Canonic RM expansions: 

The eigenvalues of the RM matrices, for the non-positive 
canonic RM expansion may be obtained as the direct product of 
the eigenvalues of the core matrices. The characteristic 
equation of the core matrix for a complemented variables being 
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an irreducible polynomial of degree 2, the roots are found 
in the extension field only and the roots occur as a set 
viz* y — a. and ■ =l+a« The number of eigenvalues of the RM 
matrices is dependent on the polarity function i.e, on the 
number of variables of the function and the number of 
variables complemented. Hence, the number of eigen-values 
of two RM matrices may be the same but the pattern of the 
polarity function being different, their eigenvectors would 
be different. 

For the polarity functions in which there is only one 
variable complemented, all roots of the characteristic 
equations for the RM matrices are to be found in the extension 
field alone. This is unique only to the case mentioned, 
because irrespective of the number of variables, if there is 
only one variable complemented of these, then the eigen— values 
of the core matrix of a complemented variable manifests 
the entire operations of the direct product and therefore 
the roots are found only in the extension field. 

In general it may be written as 

R = [J?] ® ... [i?] times 

Spectrum of R = 

[Jj ® . . . [J] © ClJaJ ® [i] ■» ti“°® 
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1(2““^)' 


m 


a 

1+a 




l+a(2”‘’’^) 


i.e. Y 
Y 


a, * 2 


,m-l 


1-Hx , * 2 


>ni~l 


Cases where more than one variable occurs complemented 
the eigenvalues y=1, y=o and Y=l+« e« obtained and satisfy 
the condition 


">1 + ">0 + “l+o “ ^ 


.m 


C3 .5.4) 


as utmost 2™ eigenvalues may be had for the RM matrix of order 
0,0*2"". AS a result of the polarity function, the first three 
variables in this category may occur 

i\} ^1^^ “ 


1(4 

a(2) ^ 
l+a(2) 


or 


hV ® ® ti3 


1+a- 


or 


[i:„] ® [i3 ® 


l+a(2) 


1(4) 

a(2) 

l+a(2) 


is being indicated in parenthesis) 


(for ease, the multiplicity 


49 


Now the variables followinoj may be true or complemented, and 
would occur as per the polarity functions, if the variable is 
true then 


a(2, 
l-i-a(2) 



1 ( 8 ) 

a(4) 

l+a(4) 


(3.5,5) 


i.e 




"i = and 


(where the prime on the multiplicity indicates the multiplicity 
after the direct product operation) and if the variable is 
complemented, then 


1(4 
a(2 
l+a(2) 




1(4) 

a(6) 

l+a(6) 


(3.5.6) 


1 . G . 


■"1 = 


ala' 


m 


1+a 




It was observed that the eigenvalues, though a function 
of the polarity function, depended more on the number of 
variables complement edj the odd number of complemented variables 
demonstrating a common feature amongst them and the even 
number of complemented variables demonstrating another 
common feature amongst them. 

Starting with two complemented variables, the resulting 
eigenvalues would be 
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C 




a . 
l+a-' 


1 ( 2 ) 

a(l 

l+a(l2 


(3.5.7) 


then with other complemented variables following this, the 

rule given in (3.4,6 l ) could be applied. It was observed 

that when even number of variables are complemented 

mi = when odd number of variables are 

complemented then “ '®l+a‘'^* which results as an 

outcome of the operation over the irreducible polynomial 
2 

Y +Y'**1*0» 1*^® eigenvalues of ten complemented variables 

are considered before generalisation. 


The eigenvalues of t complemented variables, obtained 
as a direct product of the eigenvalues of the core matrices 
for complemented variables, is tabulated to illustrate the 
observation that the eigenvalues are indeed a factor of t being 
odd or even and that the rules 

"'l ^ ^ even) 

and “ '^l+l”’^ (for t odd) holds. 

The table shows the multiplicity of each root. 
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Table 3,5,1 


No. of 
varia- 
bles 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


The 

eigen- 

values 

d 

‘ 1+a 

a 

1+a 

a 

1+a 

a 

1+a 

a 

1+0 

a 

1+a 

a 

1+a 

a 

1+a 

a 

1+a 

a 

1+a 


1 

■ 

1 

2 1 

i 

6 

m 

HE 

mm 

22 

42 

86 

170 

I 342 


a 

■ 


3 

5 

u 

21 

43 

85 

171 

S 341 


' l+a» 

■ 

B 

3 

5 

iJ 

21 

43 

85 

171 

341 



Now generalising (3,5,6), two sets of forrrulae on^, each for 
t odd and t even may be written as 


1 ) t odd: 


Y = 


Y = 


, 2^-2 
1 ^1 = 

2^'*‘^+2 

a m„ = 


a 


2t+1^2 


Y = 1+a 


(3,5.8) 


ii ) t oven; 


2^+2 


Y = 1 ^ 

Y = a 


2t+1^2 


Y = 1+a m 


2t+l_2 


1+a 


^4 




(3.5.9) 


876Dii 
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If t = m, i.e. 
the condition (3.4,4) 

_ 2K2 ^ 


all variables are complemented then 
must also hold i.e. 

2^'*'^+2 . 2"^ ■^^+2 „t ^m 

6 ^ ~ ^ 

2 !!^ + = 2^ = z™ 

O D 


Now if these t variables were complemented out of 
the total m variables of the function, the remaining (ra-t) 
variables being true, the eigenvalues of the RM matrices is 
obtained by multiplying the eigenvalues obtained for t 
variables by a factor of Eqn. (3,5.3) refers. 

Therefore the multiplicity of the eigenvalues of RM matrices 
for an m variable function, with t variables complemented as a 
function of m and t are obtained as 


i) t odd ; 


If 


, o(ni-t) 2^-2 


(3.5.10) 


Y = a = 2' —3 


- oCra-t) 

Y « 1+a = 2 


ii) t even t 

„ _ a*"'-!) 2!i2 

”1 " ^ 3 


Y = i 


(3.5.11) 
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Y = a m^=2 


(m-t) 2"*^ '^^+2 


Y = l+oc = 2 


_ o ( m - t ) 2^^^+2 


and these sets of formulae also satisfy the condition (3.5.4) 


mT+m^ + m, . = 2 ’^, i.e. 

1 a 1+a 


t+1, ^ ot+1, 


e o(ra-t) ( 2 ^ . 2 ^ ^•"+2 . _ 

.e.,z \.o f ^ T (■ ) ~ 


-.m 


2(n-^t) (2^ ^ 2^ ^ . 2 m, 


Illustrative Example 3.5.1 

The number of variable m = 3 
with polarity x^ ^ 2 X 1 ^ 

The corresponding Reed Muller matrix R = f-ii^ CH^ 


00 

10 

00 

00 

00 

11 

00 

00 

10 

10 

00 

00 

11 

11 

00 

00 

00 

10 

00 

10 

00 

11 

CO 

11 

10 

10 

10 

10 

11 

11 

11 

11 
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the eigenvalues are obtained as a direct product of 

jL 1+a ® 1 spectrum of R is (a ,a, l+a,l+a,a,a, 

l+a, 1+a ), 


For the same number of variables if the polarity now is 
made as 


the Reed Muller matrix R = 0 


00 

00 

10 

00 

00 

00 

11 

00 

00 

00 

10 

10 

00 

00 

11 

11 

10 

00 

10 

00 

11 

00 

11 

00 

10 

10 

10 

10 

11 

11 

11 

11 


The eigenvalues are obtained as 2+a^ 1^1 ^ direct 

product and the spectrum of R is (a, a, a, a, l+a,l+a,l+a,l+a) , 


It is thus seen that when there is only one variable 
being complemented, then all roots are in the extension field 
only and irrespective of the pattern of the polarity 
function, the eigenvalues are the same. 

Illustrative Example 3,5.2t 

Consider a 4 variable switching function 
with polarity x^x^x^^^, 
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the corresponding Reed-Muller matrix R = 
a direct product 


The eigenvalues are obtained as ( 

the spectrum of R being (a,a,a,a,a,oi;,a,a, l+a,l+a,H-a, l+a, 
1+a, 1+a, 1+a, 1+a) and for the same number of variables, with 
the polarities now assigned as 


The RM matrix becomes R = C^l^ ® ^11^ ® ^11^® ^11^ 


The eigenvalues are obtained as a direct product of 


The spectrum of R being (a, a, a, a, l+a, l+a, l+a, l+a, a, a, a, a, 
l+a , l+a , l+a , l+a ) . 

It is once again observed that when there is only one 
variable being complemented, then all roots are in the 
extension field only and irrespective of the pattern of the 
polarity function, the eigenvalues are the same. 

Hence fox 

m * 3, y=at and Y=l+a, and 

m « 4, Y=a, ra^=8, and y=l+a, 


I 
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Illu$trative Example 3.5,3: 

A 7 variable function with five variables complemented 
with the polarity function 1000 100 is considered. 

The function is 

( Xyj x^, x^i# and 

the eigenvalues are obtained as direct product of the eigen- 
values of the core matrices as 

Ci3® hi.] 


the number of 
eigenvalues 
after each 
direct-product 
operation 


1(2) 

a(2) 

1(4) 

1(4) 

1(8) 

1(24) 

1(40) 


0^(2) 

o((9-) 

a(6) 

a(12) 

a(20) 

a(44) 



(ft SI-) i 

a^(6‘ ) j 

a^(12) 1 

1 2^(20) 1 

10^(44) 

! ) 


and by formulae, t being odd 

for Y = nij: = 4,10 = 40 

for Y=« . \ = 4.11 = 44 

for y = !«, = 4.11 = 44. 

Therefore mj_=40, m^=A4 and mj^^^ = 44 and E = 2™. 

The same were obtained as direct product of eigenvalues 
of the core matrices also. 
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If for the same function the polarity is 1010100 
then m=7, t=4, 

for Y=1 nij^ ~ = 8.6 = 48 

for Y=a = 8.5 = 40 

for Y~l+a ^ = 8.5 = 40 

and £ = 2"^ 

3.6 CHARACTERISTIC EQUATION: 

The condition for a non-null solution to the equation 

(3.5.1) 

Ru = Y^ 

for R, a matrix and u, a vector and y> ^ scalar is that 
R - Yl is singular i.e. |R - yI! = 0. This condition is 
called as the characteristic equation of the matrix R. For 
the matrix R of order - n, the equation is a polynomial 
equation in y of degree n. The roots of this characteristic 
equations are called as the characteristic roots or the 
eigenvalue ( Secs. 2,6. and 3,5 refer). 

For the RM matrices, which are obtained as a direct- 
product of core matrices, the eigenvalues have also been 
obtained as direct product of the eigenvalues of the core 
matrices* Thejse‘’‘very eigenvalues being defined by the 
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characteristic equation, it is feasible to’ construct the 

characteristic equation ficormthe eigenvalues, if Y 2 * • 

are the different eigenvalues with multiplicity • • * 111 ^ » 

Yi Y2 Yi^ 

then the characteristic equation may be written as 

(Y“Yi) ^ (Y“Y2) ^ (Y-Yj^) ^ = 0 (3. ‘6.1) 

And in the case of RM matrices, whose roots are over the 
extension field GF(2^), with the irreducible polynomial 
Y^^Y+I » 0, and all operation being modulo 2, the equation 
may be written as 

(Y+1) ^ (y^+Y+ 1) ^ = (Y+I) ^ (Y +Yti) = 0 

(3.6.2) 

and these equations may ■'be specifically written for the 
m variables, of which t are complemented as 

(Y+1) ^ = 0 for t = 0 m^ = 2 
(Y^+Y+1) for t = i m^ = 2”^"^ 

m 1 o 

and iy+l) ^(y'^+Y+I) ^ ^ 1 < t < m 

and +/- accordingly as t is odd/even respectively. 

3,7 EIGEN" -VECTORS OF THE REED-MULLER MATRICES; 

The roots of the characteristic equations of the 
Reed-Muller matrices, being limited to the three kinds of 
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roots that may be had and hence their occurance with 
multiplicity greater than 1, suggests that it may not be 
feasible to uniquely associate an eigenvector to each root. 
Eigen-vectors corresponding to multiple roots may be found 
using the Jordan canonical form of matrices [l3]. A more 
specific method of calculating the eigen-vectors, when the 
characteristic equation has multiple roots is by the use of 
the generalised inverse of the characteristic inmatrix of 
the RM matrix R. 

The condition for which non-null solutions of the 
vector u, may be obtained for Ru = Yu is that (R-yl) is 
singular. For singular matrices, only the first of the four 
Moore-Penrose conditions for the unique inverse may be 
satisfied and hence only the generalised inverse may be 
obtained. 

Using this generalised inverses of the characteristic 
matrix and an arbitrary matrix z, the non— null solution of 
u may be obtained as 

u = [(R0 yI)" (R<©Yl)0l] z (3.7.1) 

Further, the number of such solutions that (R-yI)u=0 may 
have are determined by the order and rank of the chara- 
cteristics matrix for that root, which may be written as 

n-r (R(£)yI) (3.7.2) 
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i.e. these may linearly independent eigenvectors may be 
obtained and with this as the basis, the entire span of 
the eigen-vectors may be had. 


The method of obtaining the modal matrix, i.e, the 
set of all non-null u as the columns, the algorithm given in 
Sec* 2,6 may be used. The characteristic matrix (R-yl) may be 
partitioned as r- 


RSlTkl =Bk = 




(3.7,3) 


for each eigenvalues Yjj» where Rj^ is non-singular with the 
same rank as the characteristic matrix (R-Yj^I). And, if Rj^ 
does not satisfy this criteria, then by operating with the 
permutation matrices on the characteristic matrix, Rj^ may 
be made to fullfil this condition (Sec. 2.8 refers). Then 
the generalised inverse obtained for the characteristic.' 
matrix is 




R 


0 


-1 


0 

0 


(3.7.4) 


Further, considering the arbitrary vector Z and partitioning 
it as 

/ Z* = [-V, r'Wj' (3.7,5) 

and substituting in (3.7.1), a number of non-null solution 
governed by (3.7.2) for arbitrary w of order n~r(R-Yj^l), 


is obtained as 



61 




+Rk Ck VO 


w 


(3.7.6) 


The modal matrix so obtained, may be modified to fit 
the field over which all operations are being implemented 
to enable obtain the same with ease. The arbitrary matrix 
w in the modal matrix (3.7.6) may be qualified in an unique 
manner to ensure that the columns are indeed unique thus 
giving the n-r(H-Ykl) unique vectors, being the set of 
linearly independent eigen-vectors. The property of the 
arbitrary matrix w, therefore woulct have to be full-column 
ranked* Such a matrix in GF(2) would be an identity matrix. 
And, hence the model matrix for the Reed-iviuller matrices 
may be written as 



(R^-h(C^)(I) 



“k = 

I 

■ 

— 

I 



i 



The order of the identity matrix being n-r(R-Ykf) where 

R, ” ' is the non-singular submatrix of the characteristic 
K 

matrix of rank, r(R-Ykf) matrix Ck» another 

submatrix of the characteristic matrix whose rows are the 
rows of the Rk matrix extended and the columns are the 
columns remaining from the characteristic matrix after 
the matrix Rk has been obtained. 
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This method of obtaininc, the modal matrix may be 
applied to both simple and multiple roots and also where 
the roots are in the extension fields. 

Illustrative Example 3,7.1: 


R = 


For a 2 variables Boolean function with polarity function 11 
the FiM matrix for the positive canonic txU exoansion is 

10 0 0 
110 0 
10 10 
[1111 

and the characteristic matrix is 

0 0 0 0 
10 0 0 
10 0 0 
1110 


(h-Yl) 


J 


the order and rank of the characteristic matrix are 4 and 2 
respectively* 


B. 


R, 


D 


'k 



’*10 00 


11 10 


00 00 


00 00 
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and Rj^ being involutory I = Rj^ 

rlO] rOO j _ rO0 1 

’-ii-' 4o-^ “ ^lo-* 


and 


and matrix w, which is arbitrary of order n-r(i - yl) could 
be choosen as the identity matrix of this order, viz. 

- 10 - 


W = [-] 


r 


— 

Rk“^Cj^ w 


rOO-] rlOn 
LiqJ LqI-^ 

W 


10 



01 



— ,„n.n- 


and thus 

0 0 
1 0 
1 0 
0 1 

illustrative Example 3,7.2: 

Considering the 2 variable Boolean*” function with 
polarity function 10, The RM matrix is 


R = [i?] ® [?h = 


0100 

1100 

0101 

1111 


The characteristic matrix 
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(R-YI) 


Y 1 0 0 

1 Y^ 0 0 

0 1 Y 1 

2 

1 1 1 Y 


the rank of the characteristic matrix is 

r(R-Yl) - ^ 

Therefore there are n-r(R-Yl) LIN EVS which is 4-3-1. 
Partitioning the characteristic matrix, after permuting 
we obtain 


B = 


1 Y 
0 1 
1 1 
Y 1 


0 0 
Y 1 
1 Y^ 
0 0 


for the root y = 


r 


R, 


a a a 

1 c? 1 


(? 1 


a 


and Rj^ = 


O 

I 


0 

0 1 Y 

and Cj^ = 

1 



2 

111^ 

1 

J 

Y 


with w = [l] 


Y==a 


and hence the model matrix for y " ® 


2 “1 


0 

a a a q 



1 1 L 


0 

2 


2 

1 


a 

a i a 



(1) 

11 

1 

j 


^ — 


L. 


which indeed is an 
eigenvector. 
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for Y = (l+a) the eigenvectc?r-imay be found as 
— — — 

0 

0 

a 

1 

3.7.1 Eigen-vectors of RM Matrix, for the Positive-Canonic 
RM Expansion: 

However, another algorithm for finding the set of LIN 
eigenvector of the ElM matrix which represent the positive 
canonic Reed-Muller expansion has been developed. This 
algorithm is based on the concept of null-space of a matrix. 

From the definition of eigenvalues and eigenvectors 
given in Section 2.6, viz. 

Ru = yu for u a vector and y scalar 
it may be written as 

(R - yI)u = 0 

and all vectors u that satisfy this condition forms the 
modal matrices, with the vectors being the columns of this 
matrix. The null space of a matrix may be defined as the 
set of all column matrices x which satisfies the homogenous 
linear equations 13x = 0. An easier approach to obtain the 
null soace of the matrix is by making the matrix B full 
row— ranked and finding a matrix C which satisfies B C — I 
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and taking those columns that satisfy Bx=0, the null-space 
may be obtained. For a matrix B = (R-yl) this may be written 

as 

(3.7.8) 




such that 


B^C^©B2C2 = I 

or equivalently the same equation may be written as, 


®11 ®12 


^11 *^12 


~I 0 : 

®21 ®2^ 


Si 

j 

0 I 


(3.7.10) 


or those C^2 *^22 satisfy the conditions 

_ O P rs O (3*7#H) 

^11*^12® ®12 ^22 

®21**^12® ®22 ^22 “ ^ 


Q the 

such that C 2 = j L 

of tha matrix rt and being the null-space of the characterl- 

Stic matrix (R-yl). 

It is observed that for '=22 = 

appropriate order, the modal matrix may be formed. The 
Choice of W and I also depends on the choice of and 

which satisfy the conditions 


(3.7.12) 
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®21*^12 0®22 • ^22 “ ^ (3.7,13) 

and 

®21*^ll©®22*Sl ~ ^ 

The matrix B 2 may be chosen arbitrarily, with the condition 
that the matrix B’ is full-row ranked , Formats of matrix B, 
always suggests that the LIN vector of may be obtained 
^'^ 21 '^ 22 ^ choosing 82 ^ the null matrix and B 22 
the identity matrix of appropriate order, i.e. 


0,*Q(^ 1. 1 = I is satisfied 

and hence the resulting modal matrix takes the form M = [j], 
where Q is of the order n/2, n being the order of the Heed- 
Muller matrix R and Q also takes the same form as the^matrix 
R, and I is of the order n/2. The identity that 

Rank + Nullity = Order is also verified as a result 
of this. 


The matrix U in the modal matrix M, may also be 

defined as the matrix Q, the characteristic matrix, with 

Y==l, and is of order n /2 or the same may be obtained as the 

characteristic matrix, with Y”l» matrix fi’, obtained as 

(m- 1 ) 

a direct-product of the core matrix times and I is the 

identity matrix of order n/2. Alternatively, Q may be written 
[H* • I], The span of the column vectors of the modal-matrix 



I 
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as the basis v -ctors, the complete set of eigen-vectors of 
the RM matrix may be found. 


Illustrative Example 3.7*3; 

Given that the no. of variables in the switching 
function, m“2. 


The Cj^s for the expression 


to® Cj^tj^©c^t^0C3t3 

where t's are the minterms and c^ ^0,1^ indicates whether 

JL 

the function value is zero or one for the related input 
and as the mintermsare disjoint, the inclusive OR, 

operation is replaced by the exclusive OR operation. For 
this example let the c^s be oi 01. The other sets of c^s 
of all possible 2'® values that this may take are the 
sixteen combinations of Is and Os, from 0000 through 1111. 


R = 


The transform matrix, R for this example is 



”10 00 ” 


”10 00 

11 00 


00 10 

10 10 

or 

11 00 

11 11 


00 11 ^ 


and the eigenvalues are 

1 i.e. Y“l. ">1 
1 


% 


= 4 
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To find the modal matrix, 
of characteristic matrix, 



we like to find the null space 
which in the simplified form is 


[° and I = rl 0i 




and 


M = 


0 0 
1 0 
1 0 
0 1 


and the set of all eigen-vectors 

^0 0 0 0 
0 10 1 
0 10 1 
0 0 11 


H = 


The characteristic matrix of the transform matrix R, 

0 0 0 0 


is (R-yD = 


10 0 0 
10 0 0 
1110 


and then 


find the null-space the matrix B so obtained 


i.e. B = [f ? ? Sj and obtain the set of x. such that 
1 1 i u 


Bx s= 0 
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and B * 


000 “ 

► 

^1 **’ 

r 

0 

1 1 


^2 


0 


^3 


0 


" 4J 


0 


B. 


B , 





1000 

may 

be obtained 

as 

1110 

0010 




0001 

0 

0 


we may write 


’*‘1000* 

1110 


^11 ®12 


10 0 0 

0 10 0 

0 0 10 

0 0 0 1 

mtmt 

M 

^21 ^22 

•Ml. 

0 0 10 

0 0 0 1 






'12 ‘•11 


rlOn . ^ rOO] ^ rOOj 
Ln J ^22 ‘-lO-' *■00’' 


and c ^2 ^22 ^01^ ” ^01^ 


^22 ‘• 01 -' 


and hence Cj^2^ai^ ^01^ ^10^ “ ^ 

rlOi rlOi r00-| , rlO-i rOOn ^ rOOi 
^12 * ^ll** ’•01*' ‘•10-' '■11'' ‘■10*' '•10'' 

(as the matrix is involutary 


^2 ’ 


00 

10 

10 

01 
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Illustrative example 3,7.4: 

the 

The number of variable in / positive canonic fiM 
expansion is given as m=5, 

a) The order of RM matrix is n = 2’® = 32 

b) The modal matrix may be easily constructed as 


0000 0000 0000 0000 
1000 0000 0000 0000 
1000 0000 0000 0000 
1110 0000 0000 0000 


M=[?]= 


1000 0000 0000 0000 
1100 1000 0000 0000 
1010 1000 0000 0000 
1111 1110 0000 0000 

1000 0000 0000 0000 
1100 0000 UOOO 0000 
1010 0000 ®000 0000 
1111 0000 1110 0000 

1000 1000 1000 0000 i 
1100 1100 1100 1000 
1010 1010 1010 1000 
1111 1111 1111 1110 

1000 0000 0000 0000 
0100 0000 0000 0000 
0010 0000 0000 0000 
0001 0000 0000 0000 

0000 1000 0000 0000 
0000 0100 0000 0000 
0000 0010 0000 0000 
QOGQ 0001 0000 0000 

0000 0000 1000 0000 
0000 0000 0100 0000 
0000 0000 0010 0000 
0000 0000 0001 0000 

0000 0000 0000 1000 
0000 0000 0000 0100 
0000 0000 0000 0010 
0000 0000 0000 0001 



CHAPTER 4 


GBiilERATION OP PBBD-MJilER EXPANSION OOBEEIOIBNTS 

Developing algorithii's , which synthesise any arbitray 
switching function, using a aiinimum number of logic circuit 
has been a perpetual problem in switching theory-* Boolean 
realisation use AND/OR/NAND and NOR gates, whereas these 
functions may also be implemented using Exclusive OR gates* 

The advantages of designing with Exclusive OR gates, were 
outlined in Sec. 1.1* Applications using Bx-OR gate reali- 
sation has suffered because there are not mapy algorithms 
that may be easily applied to find the minimum solution, as 
are available in the case of Boolena vertex realisations. 

In the eiclusive-QR realisation for the minimum design, it 
is desirable to mAnimise the number of terms, in exactly 
the same manner that it is desirable to minimise the number 
Of prince complicants (ANN ) terms in a conventional Boolean. 

realisation* 

Bl«en structure of any matrlr characterises the 
matrix and therefore it is possible to aeaooiate theel«en- 
atruoture of the Heed-Mu Her matrices, obtained in the previous 
chapter, in the generation of the positive oanonlo- and non- 
positive oanonie Reed-MiHer expansion ooef f ioie nts . Thus 
the eigenstruoture of the expansions 1 * able to relate the 
mint^ ooeffloiehts (obtained fremthe troth table.- of the 
Boolean function) to positive oanonlo and non.-posltive 
fDanonle RM expansion coefficients. 
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The algorithms developed h<=re are for the two different 
rM matrices* In the algorithm, for the positive canonic Rl’i 

jK 

expansion, eigenvectors pertij^ns the complete set of minterm 
coefficients vector space C, into cosets and hy associating 
the coset leader with an eigen vector, the positive canonic 
RM expansion coefficients may he obtained* The other algorithm, 
makes use of the eigenvalue of the non~positive canonic RM 
matrices and is able to associate the minterm coefficients to 
its corresponding non-positive canonic RM expansion coeffi- 
cient for a given pt^lsr%’f unction of the Boolean function* 

Both the algorithms may be programnted, when the matrices 
cannot be handeld i.e. for Boolean functions with large number 
of variables, and may be enployed profitably in the minimisa- 
tion problems* 

4-1 BJRPOSB OP aSHBRATION OB THE NON-ROSITIVB CAITONIC RM 
EXPANSION CQBFPIOIBNTS 

The generalised Reed-mller (&RM) expansion may be 
expressed as (1*7) 

f(3c^, 3t2' » V = ^ ^2 ^ 

^ **** 

where each 5^ input variable may appear in positive or negative 
form but not in both forms- 



74 


The number of zero~value«i coefficients in the 

GrRM expansions for any given function f(x) and hence the number 
of remaining Ai^D and Exclusive terJiB in the r<»sulting reali- 
sation, will vary depending upon the chosen polarity of each 
Xj^ input variable* The extreme polarities of positive and 
negative do not contain more serovalued terms than some other 
non-positive canonic RM expansions. Hence the criterion for 
minimisation is the selection of that non— positive canonic 
expansions which has the maximum number of zero (minimm 
number of non-zero) coefficients in the expansion. A com- 
parison of all possible non-positive ogsnt>nic RM expansions are 
involved in order to determine a miniamm non-positive canonic 
RM expansion [21], [22], [23], [27]. It is difficult to 
generate these 2^ positive and non-positive canonic HM ex- 
pansions and these are not intutive either. Constraints 
on the variable being complemented may give rise to minimum 
solution, A method has been developed [7], where geometric 
plotting and grouping have been outlined and thus after manu- 
pulatiftg the coefficient values, a minimised non-positive 
canonic RM expansion may be achieved without involving 
exhasutive search procedure. A computerised in place algo- 
rithm has also been devised [IO]. The algorithm enables 
sequentially generate all non-positive canonic RM expansion 
coefficients, and the algorithm being fast, an exhaustive 
search for the minimum is feasible. 
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Ihe table below dearly shows the advantage of the 
non-positive canonic RM expansions. 


TABI® 4.1.1 


POSITIVE AMD HOU-POSITIVB 
OANOiVIC Rii EXPANSION Ii'!iPlEi''iBi'ITATIOE 


s. 

Number 

Minterm* 

Positive 

Non 

-Positive 



Polarity* 

ITo. 

of vari- 

Coeffg « 

Canonic 

RM 

Canonic 

Rm 

Pattern 


ables in 
the 

f<x^ »^2**** \ 

JjBxpansxon 

Expans non 

of the in- 
put vari- 


function 


X.0R 

ANI| : 

BOT 

^XOR 

AND 

NOT 

able 

1 

.-4 

.2528 

6 

5 

0 

2 

2 

2 

6 

2 

4 

21845 

2 

1 

0 

1 

0 

1 

1 



(101010...) 








5 

4 

65280 

1 

0 

0 

1 

0 

0 

0 



_(-Q--«P0l ...11; 








4 

4 

43690 



1 

0 

0 

1 

0 

0 

0 



(0101 ...) 









4 

255 

2 

1 

0 

1 

0_ 

1 

8 

6 

5 

1 0000 1 

16 

14 

0 

6 

B 

bb 

f 1T 



1 




5 

■eW 


i 23 


5 

20000 

i2„ 



.12_ 

-0_. 

-6. 

B 

B 

i 16 

8 

5 

1^21655765 


Li _ 

0 


0 

1 

i .1 

mm 

5 

227216640 

16., 

HH 


7 

U.6.. 

.2 

: .14 

[q_ 

5 

1671935 

,-2^, 

1 

.Oj 

B 

0 

J 

8 

9 

6 

1 1 0000000000 

2P.. 




Bb 


54 

!|^ 

6 

^1241814220 

24- 

.20.-. 

0 

12 

12 

^4—. 

60 


6 

l2212I624128 

26- 

.22.. 

0 

14. 

-14— 


61 

9B 

6 

.1224507890 

20„. 

18 

0 

10 

10 

1 

22 

L5_ 

6 

.102039405,0 

28.. 

.24._, 

.0.. 

.11_ 

-11— 

_2 

40 


* Decimal value of the binaryrepresentation of the minterm co- 
effiojants and the polarity function. When the pattern of min- 
term coefficients is not indicated, then it is a random pattern. 
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Minimisation criterion, depending on the realisa- 
tion may "be applied ; one is to minimise th» number of ACfD 
gates, another one is to minimise the overall number of inputs 
to the Ex— OR gates or to minimise the overall number of 
inputs. However, the choice of criterion being the economical 
imp 1'- mentation of the situation to which the d'^sign is to b^^ 
applied. The penalty to be paid for the non-positive canonic 
RM expansion implementation is th'^ sub-optimum number of 
product terms, while the numb'^r of input variabl-^-s or the 
number of inverters required is gen-rally less than for the 
minimum form [jO]. 

4,2 G-EiTSHATIOH OP POSITIVE CiJNOi'JlC RM EXPANSION OOEPE ICIBHTS 
USING BIGBiWAlUBS 

The vector*' space of all combinations of possible 
values that c^ s may take get partioned by the eigenvectors of 
the transform into subgroups, and the RM coefficients for a 
particular subgroup are governed by a particular eigenvector. 
Thus it IS observed that this vector space is partioned into 
cosets and the elements of the subgroup or the cos^t is in- 
variant i.e. the RM expansion coefficient is found from within 
the coset, govern<='d by an eigenvector. 

This if h^, hg, \ are the eigenvectors of the 

transform, then cosets may be formed as shown (all op<=>rations 
are modulo 2 addition) 
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I 2 . 

®2 0^1 ® 2©^2 ^ 20^3 '■*’ 



® k ©^1 ^ k ©^2 ' S ] i :®^3 •*•* ® k ®\ 

and in this all ® lecome to coset leader. While 

forming the oosets, instead of g^0h^ being any element that 
has not appeared earlier, it is the n<=xt low-^st element after 
h^ , that has not appeared earlier. Thus the coset leader is 
the minimum vector in the coset. In the algorithm hit-hy-hit 

modulo 2 addition of all eigenvecboi^ to the c^s gives the 

\ 

entire coset and by chocalng the minimum vector the coset* 
leader is thus obtained. 


How after identifying tht coset leader, the eigen- 
vector that relates the c^s to the M expansion coefficient 
may he found by multiplying the coset leader with the charac- 
teristic matrix of the transform. Having found the eigen- 
vector that relates the o^s to the expansion coefficient, 
a mere bit-by-bit adaition of c^ to the eigenvector gives the 
RM expansion coeff icaents . 



i 

Qi • 




OJ 

» 


•H 

pH 


The algor it ha\ is schematically depicted above. 
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Illustrative Example 4,2*1 

The cosets for m = 2 could he written as shown helow 
(the decimal value of Cj^s are given for ease of notation). 


1 Ooset b 

i Leader*^ 

3 

2 

1 




0 

6 

8 

14 

Eigenvectors 

4- 

1 

1 

7 

9 

15 

14 

2 

2 

4 

10 

12 

8 

3 i 

3 

5 

1 1 

13 

i 

6 


Thus for the given c^ 0101, which is notationaly 10, 
(Ten) hy car' ying out the hit-hy-hit operation on all eigen 
vectors we get 10, 12, 2 and 4, the minimum of which is 0100, 
which is 2, the coset leader also, Now carrying out a vector 
matrix multiplication, i.e. the coset leader and the charact- 
eristic matrix. 


“o' 


"o 

0 

0 



'o ! 

t 

0 

— 

1 

0 

0 

0 


1 '' 

0 

j 

1 

0 

0 

0 


0 

I 

1 i 


1 

1 

1 

0 ' 


0 

_ _ 


we obtain the RM expansion coefficient determining eigen 
vector, which is 0001, and bit-by-bit addition of this to the 
c^s gives the RM expansion coefficient. 
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i.e. 


r 

f 0 
! 0 


! 




+ 



1 i 

! wh-xcli ia so. 


Illustrative Example 4,2,2 

The coset formation for a 3 variable function is shown 
in the table. The coset formation contains all the 2^ 
possible combination which the set o-f o^s , the minterm coeffi- 
cients could take. Hence once the table is formulated and 
the appropriate coset leader associated with th® corresponding 
eigenvector, using the algorithmj the table is available as 
a ready reference or as a mere look-up table. The effort 
therefore is only to do a bit-by-bit addition, modulo 2, of 
the minterm coefficients with the corresponding eigenvector 
to give the Reed-Muller expansion coefficients. However, 
in this table the positive canonic RM expansion coefficients, 
corresponding to each minterm coefficient is also given as 
the denominator to it i.e, after the vector obtained as a 
result of bit-by-bit addition of the minterm coefficient 
with the positive RM expansion coefficient determining eigen- 
vector found from the table. 



TABLE 4.2.1 

C OSE'T LEAEE E, BiaEWECTOR LOOK~aP TAB LE 



The above table also gives the positive cationic RH expansion coefficiints 
The denominator O'f each of the ttiinterm coefficient entfy is the posijive 
canonic RM expansion coefficient corresponding to it. 
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Illustrative Example 4,2*3 

The number of variables in the positive canonic Ei^x 
expansion is given as m =* 4 , 

(a) The order of the RM matrix is n = 2 ® = 16 

(b) The modal matrix may be obtained as 


E = [«] = 


00000000 
1 0000000 
1 0000000 I 
1 1 1 00000 
1 0000000 
1 1 001 000 
10101000 i 
11111110 1 
1 0000000 . 
01 000000 I 
001 00000 ! 
00010000 
00001 000 
00000100 
oooiool 0 
00000001 


The above linearly iadep-=-adent eigenvectors as the 
basis, the span of these vectors may be found to be consis- 
ting of the 256 eigenvectors of the RM matrix. There are 2^ 

1 6 

possible combinations of the minterm coefficienta i.e. 2 =6553^^ 

So Imrtead of handling this many combinations, we may store only 



the 256 eigenvectors of the RM aiatrix and hy hit-by-bit addi- 
tion of the given minterai coefficient to these 256 eigen- 
vectors and picking the aiiniimii? vector, the coset leader is 
obtained. A 

The process of associating the expansion coefficient 
deteraiining eigenvector to the coset leader could be available 
as a look— up table, as this is independent of any operations 
with ttiinterffl coefficient and therefore once the cosct leader 
is obtained, then by identifying the expansion coefficient 
determining eigenvector, bit-by-bit addition of the minterni 
coefficients gives the corresponding RM expansion coefficients, 

This algorithm is very efficient, provided the two 
look-up tables are available before hand viz. the set of all 
eigenvectors and the coset-leader to eigenvector table, then 
with only two sets of bit-by-bit additions, the EH expansion 
coefficient may be obtained. In general, 2^^^+ 1 bit-by-bit 
addition and the process of picking the coset leader is ade- 
quate and the cosot leaders could also be stored in memory 
for comparision and thus speed-up the operations. The number 
of coset leader being the total number of eigenvectors itself. 

4.3 Q-BJTBRATIOi'f OR RON-PCB ITIITB OANQtflC RH EXPANSION OQSRRICIt; 

ENTS, USIIG BIGBmXIJBS 

To enable relatflthe c^s to the positive and non-posi- 
tive canonic RH expansion an algorithm has been developed. 
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which is schematically showa below 





7 



— iv 

L 

Characteristic | 
•matrix i- 

i 

¥ y s 


{ 

! 

1 


Goeff. 
jM detera\xa 
lag . 
2.ex!jLQr_. 






GRi-i 

coeff ac i^nt 


This algorithm may be seea as a version of the oae 
presented in the earlier section and this being applicable 
to the more general case may be applied to the ’positive ’ 
polarity RM expansion coefficients also. 

This algorithm may be written as 

l%y y (R + Y I)] -[X] f [X] = RI 

There are three kind of eigenvalues that are obtainable from 
the RM matrices viz, y = 1 » Y = and y = 1 + “J also it was 
mentioned in section 3,2 that whenever y = a occurs then 
y = 1 + a also occurs. The total number of eigenvalues is 
2°^ and thds the number of eigenvalues that appear for y = 1 , 
if they occur, is always even, Bfe must be remembered that 
all arithmetic is modulo 2 and therefore the summation of 
y(R-Yl) for all ys , as far as y = 1 is concerned may be omitted 
and evaluating for y = a and y = 1 + oc, we obtain. 


IHS = [a(R + al) + (r + a )(R + (1 + «) I)][x] + [x] 
- [c<R + I+(1 + u) R+(1 +u)^ l][x] + [x] 
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= [a® + 1+ aI + R+ 0©. + al] [x] + [x] 

= [I + R] [X] + [X] 

= X + RX + X = RX = RHS QBD. 

It is observed that "weightiug the characteristic 
mat^lx-by tha root rasalts in the characteristics matrix, 
which is what has been used in the pre'vions algorithm. 

In this algorithm instead of obtaining the eigenvectors that 
determines the RM expansion coefficients, some arbitrary 
vector is obtained, which determines the GRH expansion coe- 
fficients. 


Illustrative Example 4*3* t 


Eor a 3 variah^eBoolean, function of polarity x^ x^ x^ 
i,e, m = 3, polarity = 101 


/. R 



X 


Ml'* 


X 


[lOj 
4 1 j 


fool 00000 
001 10000 
1 01 00000 
1111 0000 
001 0001 0 
001 1 001 1 
1 01 01 01 0 
11111111 


1 0000000 
001 00000 
00001 000 
0000001 0 
1 1 000000 
001 10000 
00001 W 
0000001 1 


01 000000 
0001 0000 
00000100 
00000001 
1 1 000000 
001 1 0000 
00001 100 
0000001 1 

'0m 


1 0000000 
001 00000 
00001 000 
0000001 0 
11000000 
001 1 0000 
00001 1 00 
, 0000001 1 


for the given set of say 10011001, we may obtain the ORM 


expansion coefficients as 


a = R ' 0 
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^I'bei'iia'fc© ly I w© know th.a't tliPr© is only one va.riabls being 
comp lamented i*e t = 1 and we may obbain the number of roots 
in tbe ground "field and extension "field eitber as direct 
product of the individual core matrices such as 

[]] S' [,+“j [,’] 

l.e. Y = « ; "V = 4 

Y = 1 + a j = 4 

or using the formula we have 

m = 3, t = 1 therefore there are no roots in the 

ground field and all roots are in the extent ion field and 

there are d = 2^^' of them i.e, 8, 4 of them are y = a and the 

2 

other 4 are y = 1+a (or cc ) 

The characteristic equation is 
2 

(y t Y + 1 ) “ ^ 

i.e. (y^ + Y + 1)"^ = 

Now using the algorithm the G-RM expansion coe"ff icients 
can be found as follows : 

a = [^ y(R 
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oc 

0 

1 

0 

0 

0 

0 

0 


1? 

0 

1 

0 

0 

0 0 

0 



! ~ 

’^1 

0 

a 

1 

1 

0 

0 

0 

0 


0 

a2 

1 

1 

0 

0 0 

0 




f 

G 




r\ 













2 


2 

1 

0 

a 

^0 

0 

0 

0 

0 


1 

0 

a 

0 

0 

0 0 

0 





Or, 


















3 


3 

1 

1 

1 


0 

0 

0 

0 


1 

1 

1 

a 

0 

0 0 

0 



c , 


c . 
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4 


4 

0 

0 

1 

0 

a 

0 

1 

0 
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0 

1 

0 

c? 

•0 t 

0 

1 

1 1 
1 

I 

1 

1=5 1 

“5 

0 

0 

1 

1 

0 

a 

1 

1 


0 

0 

1 

1 

0 

a^1 

1 


i 

1 

[ 

°6 ! 



1 

0 

1 

0 

1 

0 

o?Ol 

1 

1 

0 

1 

0 

1 

0 a 

0 

i 

! 

“tI 

1 

1 

1 

I 

°7 

1 

1 

1 

1 

1 

1 

1 

« 


1 

1 

1 

1 

1 

1 1 

f 



1 

1 

1 

i 

— 1 








— — « 








j 


1 


and simplifying over the irreducible polynomial. 


“~1 01 00000 “ 


"Ci~ 


■ 

01 1 1 0000 


°2 



1 0000000 


°3 


°3 

1 f1 00000 


°4 

@ 

°4 

001 01 01 0 


°5 


°5 

001 1 01 1 1 


Cg i 


°6 

1 0101000 


°7 

i 

Orj 

11111110 

1 

°8 

1 



02^0^4 


C^C^®C^Crj^8 
c ^0 ^ Y 



1 


°2 


°3 

0 

°4 


C 5 


°6 


°7 

) 

J 



winch for the givea set of c^s viz. 10011001, is 
a® = [01100000] 

and thus the required GRh expansion coefficients have bpen 
found for the polarity 101. 

4.4 GENERATION OE THB Rh EXPANSION COBEPiCiBNrS BASED ON GOOD'S 
EACTORISATION OE THE RM NIATRIGES 

An »m* variable Boolean funcbion could have any 
polarity for specified optimuffi solution. Aa already mentioned 
’in~place' computer algorithm [lO] is available for this. 

The program given at Appendix too, though not 'in-place', 
gives the same results, with the flexibility that even for any 
random polarity the non-positive canonic RN expansion coeffi- 
cients may be had in a single run. The Reed-huller matrix 
is built-up as Kronecker product and then factorised as two 
kinds of matrices using Good’s technique, which when sequen- 
tially processed gives the desired positive /non— positive canonic 
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RM expansion coefficients. As the representation in the 
program for the matrix is the existence or otherwise of an 
• 1 > element in the matrix, the computation as also fast i.e. 
only the location of the elements in the matrices determines 
the output vector. The Reod-Muller mebrix is thus factories -d 
as 




with R*. taking one of the two formats d'^pending on wheth-'r 
the variable is true or complemented. The program in its 
present form scans all the possible 2^' polarity functions in 
a binary sequence and from this the minimum number of Ai^D gates 
or inputs to Exclusive OR gates or with any oth^r specifications, 
the desired optimum polarity may be had. The limitation of 
this program is the capability of the comput«r/handle matrices. 
The program is given at Appendix G> • 


The minterm coefficients for a 5 variable Boolean 
function la given an the input. These minted coefficients 
are the same used in [10] for the erasple there. The output 
may he had in two forms, one as a sun®ary of the numher of 
gates required for each polarity function which is reproduced 


as a table hoi*®* 
gives the details 


the other besides giving the ebove data 
Of RM expansion coefficients, a^s for each 



case. The program iji its present form does not give the 
exact number of Bx—OR gates required, though the same may be 
construed and neither does it indicate the type of the 
gates required i,?, hovf many input AMD gates. However^ an 
analysis of a^^s , the RM expansion coefficients would divulge 
the number of inputs to each AMD gate, 

TABIE 4 . 4.1 


RM BXPAMSIOM COBFPICIEM'IS FOR DIBPBEBMT PODARITT OP THB 

VARlABiES OP A BOOlBAM PUMCTIOM 


B-J 
Super- 
s cript 

Number of 
input to 
BxOR gates 

Number of 
AND gates 

Numbpr of 1 

Inv'rtors 1 

(MOI gates) 

Re marks 

0 

16 

14 

0 


1 

14 

13 

1 


2 

12 

11 

1 


3 

14 

1 2 

2 


4 

9 

9 

1 


5 

1 1 

11 

2 


6 

8 

8 

2 


7 

13 

12 

3 


8 

17 

14 

1 


9 

16 

14 

2 


10 

12 

10 

2 


1 1 

16 

13 

3 


12 

9 

8 

2 


13 

10 

10 

3 


14 

7 

6 

3 


15 

11 

10 

4 


16 

20 

17 

1 


17 

18 

15 

2 


18 

18 

15 

2 i 


1 9 

16 

: 14 

3 ! 


20 

13 


2 


21 

17 

15 

3 


22 

14 

13 

3 


23 

21 

17 

4 


24 

1 6 i 

14 

2 


25 

16 

13 

3 


26 

14 

12 

5 


27 

14 

n 

4 


28 

10 

9 

3 


29 

12 

12 

4 


30 1 

10 I 

9 

4 

nr 


31 

14 i 

12 

5 
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4*5 G-SHERATI Oli Qp FilHTERjyi CO[BPPlCISi\riS EROii a G-IYEE RFi 
BXPMSIO iY CCffiPElOlBi-IIS , FOR THE GIYEF POLARITY 

The lirit-ar transformation from the inint^rai coef'fi— 
cients to the RM expansion coefficients was represented by 
the RFx Jnatrices . The expression being 

a = Rc (4.5.1) 

Therefore from the given set of non-positive canonic RFi 
expansion coefficients, the minfe^^rm coeff ici'=nts may be 
obtained as 

c = r"*^ a ( 4*5.2) 


The RM matrices were obtained as direct product of core 
matrices, and structure of the inv'rse of the RK matrices is 
also obtained as direct product of the inverses of the core 
matrices. The inverse of the core matrix is the core matrix 
itself, the core matrix being involutary. The inv'-rse of the 
core matrix, for a complemented variable is 


[0 1]*"^ = [} h (4.5.3) 

acid iuvors© of “ths RH Jiia“brioos aro olDtain d as , 

dii*pci/"*p 1 * 0 duct of tiisso iuvoi'sos coE'ibxn. 0 d as poT ttic polai*ity 
fucotiou of th6 input Thus for a polarity 

function 101001, the R ^ is 
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Illustrative Example 4,5,1 

E'er tho 3 variable Boolean function considered 
earlier, with polarity function HO, the non-positive canonic 
RM expansion coefficients were obtained as 

[a]^ = [01100000] 

Tho matrix r“^ may be obtained as 

= [ ’ °] ® [J °] [’ i] 

Therefore the minterm ooeffici-nts Q\ay be obtained as 


1 1 000000 


0 ! 

1 0000000 


1 ! 

1111 0000 


1 1 

1 01 00000 ! 


° i 

1 1 001 1 00 


0 

1 

1 0001 000 


° 1 

11111111 


0 ' 

10101010 


0 ! 
_ J 


which i$ [cj® = [10011001] 

These were the minterm coefficients originally assumed 
in the example and therefore r"^ may be obtain-H as mentioned 
in this section. 



CHAPTEf. 5 


CONCLUSIONS AND SUGGESTIONS FOR FURTHER RESEARCH 
5 . 1 3Jm JIY OF THE THESIS; 

Minimising the number of logical circuit elements 
has always been a problem. Various techniques such as 
Karnangh map and Onine McClusby's technique, and the 

recent mapping* help achieve a limited aim. The 

search was for more economical realisation and therefore 
designing only with exclusive-OR and /,ND cates were 
adopted, ( in addition, inverti2&ji. too would be required). 

This gave rise to the problem of minimisation of the Ex-OR gates 
gates and the parameters on which this minimisation has 
been defined is the polarity function of the variables of 
the [3oolean function. So the iboolean function, via its 
minterm coefficients could be related to which came ^to be 
known as the tiM expansion coefficients. The transformation 
was accordingly termed as the Flivt matrix. 

The RM matrices were obtained on direct-product of 
core matrix and thus the eigt;n-values could be easily 
obtained as direct— product . Ihe problem of multiple roots 
of the characteristic equation, ^specially, in the case oj. 
the RM matrix for the positive canonic RiV expansion was 
simplified and the eigen-values of the RM matrices for the 
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non—posi tiv6 canonic liM expansions were formulated, again 
as a function of the polarity function of the input variables. 
As the FIM matrices could be obtained as direct-product of 
core-matrices, they could also be factorised by Good's 
technique. The main advantage of this factorisation has 
been that the few non-redundant entries in a matrix could 
be made fewer, resulting in fewer number of operations for the 
transformation implementation. This aspect has been exploited 
utmost for computer implementation of the generation of 
the HM expansion coefficients. 

The purpose of generation of RM expansion coefficients 
itself was considered essential for minimisation of the 
Ex-OH designs. The eigen-structure of the RIvl matrices were 
used for the g<-neration of the RM expansion coefficients. 
Good's technique of factorisation was also put to use. 

The main dark area in the use of tx-Oh design has 
been that not many usable algorithms are available. This 
thesis is an effort to contribute towards understanding the 
problem and try and give an easier approach to it. Table 
4 . 1 .l^Chapter 4 would evince the circuit designers into 
thinking designing chips, thaugh may be exclusively 

dedicated, but extremely small in size. However, the 
economical production of the 'chips' would be the deciding 


factor 
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5.2 SUGGESTION FOR FUTHLh RESE/I^CH; 

The eigcinvectors of the Rm inacriaes, for the positive 
canonic Fif/i expansion, could be used to form lookup rabies 
and thus by a mere bit-by-bit addition, of the mintorm 
coefficients with the expansion coefficient aetermining 
eigenvector, the expansion coefficient could be obtained. 
Similar lookup tables could be formed for the FOh matrices, 
for the non-positive canonic Rivl expansions, so that the 
non-positive RM expansion coefficients too could be obtained 
in a simple manner. 

Practical designing of the circuits, using axclusive-OR 
gates, AND gates and invertors could be profitably 
worked out. The algorithm developed in this thesis or 
[lO ] , could be put to use. 

Lastly, every digital information processing essentially 
being represented by a set n iunction of m vaxiables on the. 
GF(2), such a set of functions could be reduced to one 
polynomial of one variable over the extension field, GF(2’^) . 
Such polynomials have remarkable propertlee. which may 
bo used effectively in designing of switching circuits [28]. 
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